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Simultaneous measurement of multiple qubits stored in hyperfine levels of trapped 
111 Cd + ions is realized with an intensified charge-coupled device (CCD) imager. A 
general theory of fluorescence detection for hyperfine qubits is presented and applied 
to experimental data. The use of an imager for photon detection allows for multiple 
qubit state measurement with detection fidelities of greater than 98%. Improvements in 
readout speed and fidelity are discussed in the context of scalable quantum computation 
architectures. 
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1 Introduction 

Trapped atomic ions represent a promising method for implementing universal quantum com- 
putation. Trapped ions have already met most of the DiVincenzo requirements for quantum 
computation and recent interest has focused on improvements in entangling gate fidelity 
[3E1I1IE], producing large entangled states [HI EL ancl scaling traps to larger numbers of ions 

EiEinnini]. 

In this paper we discuss the important requirement of a multi-qubit measurement capa- 
bility. In ions, state detection is accomplished by applying polarized laser light resonant with 
a cycling transition for one of the qubit states and off-resonant for the other state. The two 
states are then distinguishable as "bright" and "dark" via this state-dependent fluorescence 
|12l 1181 1141 H3j . Typical schemes collect this fluorescence using fast lenses and detect pho- 
tons using a standard photon-counting device such as a photo-multiplier tube (PMT) or an 
avalanche photo-diode (APD) . The relatively high detection efficiency of PMTs or APDs aids 
detection, but for detecting more than one ion their lack of spatial resolution means that 
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certain qubit states are indistinguishable, e.g. one bright ion out of two does not determine a 
particular ion's state. Distinguishable individual qubit state detection is particularly crucial 
for tomographic density matrix reconstruction |lfil l%|. quantum algorithms |17II18| . quantum 
error correction EH], and cluster state quantum computation [201 ED Separating the ions 
with shuttling [HJ EH EH EH or tightly focussing the detection beam EB1 can distinguish 
the qubits, but the additional time necessary for detection, possible decoherence associated 
with shuttling, and technical difficulties make these schemes less desirable for large numbers 
of ions. 

In this paper we discuss the use of an intensified charge-coupled device (CCD) as a photon- 
counting imager for simultaneously detecting multiple qubit states with high efficiency. We 
theoretically model the detection fidelity of qubits stored in S\/2 hyperfine states of alkali- like 
ions, where one of the qubit states has a closed transition to the excited electronic P state 
manifold (applicable to odd isotopes of Be + , Mg + , Zn + , Cd + , Hg + , and Yb + ). We present 
data for the detection of several m Cd + ions using a CCD imager, and discuss technical 
features and limitations of current CCD technology. We finish with a discussion of future 
improvements and prospects for integration with scalable quantum computation architectures. 

2 Detection Theory 

2.1 Basic detection method 

There are two classes of alkali-like atomic ions that are amenable to high-fidelity S\ / 2 hyperfine- 
state qubit detection. Ions that do not have a closed transition to the excited electronic P 
state require shelving of one of the hyperfine qubit states to a low-lying metastable electronic 
D state (odd isotopes of Ca + , Sr + , and Ba + ). The detection efficiencies in this case can be 
very high; typically this method requires a narrowband laser source for high-fidelity shelving 
|22l I23| although recent work using rapid adiabatic passage may relax this laser requirement 
[21] . Alternatively, one can obtain moderate detection efficiency by using coherent population 
trapping to optically shelve a particular spin state [2Sj- Ions that possess a closed transition 
to the excited electronic P state (odd isotopes of Be + , Mg + , Zn + , Cd + , Hg + , and Yb + ) can 
be detected directly, and will be the focus of this paper. Throughout this paper we assume 
that the Zeeman splitting is small compared to the hyperfine splittings. 

There are two basic schemes for this direct state detection as outlined in figureEl For both 
methods, the qubit is stored in the hyperfine levels of the Si 12 manifold with hyperfine splitting 
<^hfs- Discussing the general case first (fig-EH)? if we write the states in the \F, tuf) basis with 
/ the nuclear spin, the Si/ 2 \I + 1/2, / + 1/2) = |1) state exhibits a closed "cycling" transition 
to the P3/2I/ + 3/2, / + 3/2) state when resonant <r + -polarized laser light is applied? If the 
qubit is in the |1) state then the resonant laser light induces a large amount of fluorescence. 
When a portion of these photons are collected and counted on a photon-counting device, 
a histogram of their distribution follows a Poissonian distribution with a mean number of 
collected photons that is determined by the laser intensity and application time, the upper- 
state radiative linewidth 7, and the photon collection efficiency of the detection system. In 
contrast, when the qubit is in the S1/2 \I — 1/2, 1 — 1/2) = |0) state the laser radiation is no 
longer resonant with the transition to any excited state. The nearest allowed transition is to 
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P 3 / 2 |/ + 1/2,1 + 1/2) which is detuned by A = whfs — ^hfp, where whfp is the hyperfme 
splitting of the P3/2 states, so an ion in the |0) state scatters virtually no photons. Thus, 
we can determine the qubit's state with high fidelity by applying er + -polarized laser radiation 
resonant with the cycling transition and counting the number of photons that arrive at the 
detector. 

For ions with isotopes that have nuclear spin 1=1/2 (Cd + , Hg + , and Yb + ), there is 
another possible state-dependent fluorescence detection mechanism by coupling to the P1/2 
manifold (fig-Hs). If we apply all polarizations of laser light (cr + , 7r, and a~) resonant with 
the 5*1/2 \F = 1) — * Pi/i\F = 0) transition then the only allowed decay from the excited state 
is back to the F = 1 levels of S1/2 126] , forming a closed cycling transition. If the ion begins 
in the state |1,0) = |1), the ion will fluoresce many photons under this laser stimulation and 
we can collect these photons as above. Conversely, the state |0, 0) = |0) will scatter virtually 
no photons under this laser light because it is off-resonant from its only allowed transition to 
the P1/2 \F = 1) levels by A' = whfs + w hfp > where Whfp is the hyperfme splitting of the Pi /2 
levels. Note that to avoid an optically-pumped dark state formed by a coherent superposition 
of S1/2II, — 1), 1 1 , 0) , and |1,1) it is necessary to modulate the laser polarization or use a 
magnetic field to induce a well-chosen Zeeman splitting 27 . 

In the following sections we present a general theory of this state-dependent fluorescence by 
determining the off-resonant coupling between the qubit states. We quantify these detection 
errors in order to calculate the fidelity of qubit state detection for various photon detection 
efficiencies. 

2.2 Statistics: dark — > bright leakage 

For both the general and the 1=1/2 specific detection methods, qubits in the dark state can 
leak onto the bright transition by off-resonantly coupling to the wrong hyperfme excited level 
during detection. Rate equations describing this off-resonant pumping yield an exponential 
probability distribution of remaining in the dark state as a function of time. Once in the bright 
state, the collected photons from the closed transition obey Poissonian statistics. Therefore, 
for a qubit initially in the dark state, we expect the distribution of emitted photons to be a 
convolution of Poissonian and exponential distributions 28, 29 , as we now derive. 
The probability of leaving the dark state at a time t is given by: 

f{t)dt = —e~^dt (1) 

TLl 

where tl\ is the average leak time of the dark state onto the closed transition. Also, the 
average number of collected photons for a qubit that starts dark but is pumped to a bright 
state at time t is: 

A(*) = (l-— )A (2) 

where td is the detection time and Ao is the mean number of counted photons when starting 
in the bright state. 

We want to transform from a probability distribution f(t)dt to a probability distribution 
of Poissonian means g(X)dX so we use eqn.|2to get t(A) and then substitute into eqn.^ This 
yields the probability of the dark qubit state producing a Poissonian distribution of collected 
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Fig. 1. Relevant energy levels used in fluorescence detection of qubits stored in the hypcrfine levels 
of the S1/2 ground state of atoms with a single valence electron, with no relevant low- lying excited 
states below the excited P manifold. Energy levels are labeled by the \F, mp ) quantum numbers 
of total angular momentum and the energy splittings are not to scale, a) Detection through 
P 3/2 level with qubit stored in the S 1/2 \I + 1/2, / + 1/2) = |1) and S 1/2 \I - 1/2,1 - 1/2) = |0) 
"stretch" hypcrfine states, for any nonzero nuclear spin I. By applying <r+-polarizcd laser radiation 
resonant with the |1) — > P%/2\I + 3/2, / + 3/2) cycling transition, qubit state |1) results in strong 
fluorescence, while qubit state |0) is nearly dark owing to a detuning of A = ojhfs — <^HFP S> 7 
to the nearest resonance, where ojhfs and ojhfp are the hyperfine splittings of the S1/2 and P3/2 
states and 7 is the radiative linewidth of the P3/2 state, b) Detection through the P1/2 level with 
qubit stored in the 1 lj 0) = |1) and 51/210,0) = |0) "clock" hypcrfine states for the special 
case of nuclear spin / = 1/2. Here, applying all polarizations of laser radiation resonant with the 
|1) — ► -P1/2IO, 0) transition results in strong fluorescence, while qubit state |0) is nearly dark owing 
to a detuning of A' = <^hfs + "^hfp T' *° n sarest resonance, where ojjjfs an d ^hfp are 
the hyperfine splittings of the 51/2 an d Pl/2 states and 7' is the radiative linewidth of the P1/2 
state. 
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photons with mean A: 

ii e (x-x ) ai / vdX X>Q 



*(A)dA = ( ;i aiXo/v X — ® 

where rj = rjD^-T is the total photon collection efficiency determined by the detector ef- 
ficiency (??£>), the solid angle of collection (4^), and the optical transmission from the ion 
to the detector (T); a± = J L "x 1S * ne l ea ^ probability per emitted photon; and the A = 
discontinuity is necessary to account for the fraction that do not leave the dark state. 

Therefore, the probability of detecting n photons when starting in the dark state is the 
convolution of g(X) with the Poissonian distribution P(n\X) = e * : 

/•Ao — \\n 

Pdark(n) = 5 n e~ a ^ + / £— £-2i e (A-*oW»dA (4) 

Je n - V 

with 8 n the Kronecker delta function and e — > 0. Re- writing in terms of the incomplete 
Gamma function wc obtain: 



Pdark(n) - e- a ^^ 



(5) 



where V(a,x) = ( a \y f e~ v y a ~ 1 dy is the standard definition of the incomplete Gamma 
function normalized such that V(a,oo) = 1. 

Figure^ shows this probability distribution for various values of the leakage parameter, 
oti. Later, we will calculate physical values of ai and Ao from the atomic parameters. 

2.3 Statistics: bright — ► dark leakage 

Similarly, the bright state can be optically pumped into the dark state by off-resonant coupling 
to the wrong excited hyperfinc level. When using the -P3/2 manifold for detection this coupling 
can only occur for imperfect laser polarizations, while for the specialized / = 1/2, P1/2 case 
this off- resonant coupling is always present (since all laser polarizations are applied). As 
before, the overall photon probability distribution will be a convolution of Poissonian and 
exponential distributions, but now reversed: after some time scattering photons on the closed 
transition the qubit is pumped into the dark state and emits no more photons (neglecting 
the 2nd order effect of then re-pumping from the dark state back to the bright state). If 
we define the leak probability per emitted photon 011 = where r^2 is the average leak 

time from the cycling transition into the dark state, then with similar statistics as before, the 
probability of detecting n photons when starting in the bright state is: 

e -(l+Q2/r/)A in a /„ 

Pbn gh t(n) = - { °- + (1 + J^ )w+1 P(n + 1,(1 + a 2 /r,)Xo) (6) 

where the first term is the Poissonian distribution from never leaving the closed transition 
and the second term is the smearing of the distribution from pumping to the dark state (note 
that ati 7^ 0:2). This distribution is shown in figure[2b for varying values of cti. 

One can further separate the dark and bright distributions by using ancilla qubits and 
entangling gates before measurement |3U| . For ions with large leak probabilities a\ or 012 
this technique may be useful because it can increase qubit detection fidelity to as high as the 
entangling gate fidelity. 
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Fig. 2. Theoretical photon count histograms for ions that start in: a) the nominally "dark" |0) 
state; and b) the nominally "bright" |1) state. Both plots use Ao = 12 with varying values of: a) 
cti/r] (the leakage probability from the dark state to the bright state per detected photon) and b) 
ail r) (the leakage probability from the bright state to the dark state per detected photon). 



2-4 Calculating the atomic parameters 

2.4-1 General case using P 3 / 2 

If we assume a detuning 6 of the detection beam from the cycling transition resonance then 
for ideal cr + -polarized laser radiation the mean number of photons detected from a qubit in 
|1) is: 

A = T D T] 2 - (7) 



1 + S + 



(?) 



where td is the detection time; 77 is the total end-to-end efficiency of photon detection; 7 is 
the radiative linewidth of the P 3 / 2 state; and s = I/I sa t is the laser-ion saturation parameter. 

Next we can calculate a\, the probability per emitted photon of leakage from the dark state 
into the bright state. The dark state can off-resonantly couple to the ^3/2!^+ 1/2,/+ 1/2) 

level with probability that goes as: ~ s (^at) > w h ere Ai = cl>hfs~ ^hfp is the S1/2 hypcrfine 
splitting minus the P 3 / 2 hyperfine splitting between \I + 3/2, / + 3/2) and \I + 1/2,1 + 1/2). 
We also need to take into account the branching ratio, Mi, to determine the transition rate 
between the dark state and the cycling transition. Thus, we can write the value of a\ as: 

ai =M l( l + S+ (^) 2 )(^) 2 . (8) 



We can calculate the braching ratio, Mi, by considering off-resonant coupling between 
the dark state (| J — 1/2,/— 1/2)) and the state P 3/2 \I + 1/2, 1 + 1/2). This state can then 
decay into the bright manifold ground states |/+ l/2,/± 1/2). The generalized formula can 
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be written as: 



7+1/2 



Mx = C(I —1/2,1+ 1/2; I — 1/2, J + 1/2) V C(I + 1/2, 1 + 1/2; / + 1/2, i) 



i=/-l/2 



4/(3 + 21) 
9(1 + 2I) 2 ' 



(9) 



where we have used that C(F, F'; /, /') is the square of the Clebsch-Gordon coefficient between 
two states F — + F 1 and mp = f — ► /' given by |31| : 



C(F, F'\ /, /') = [(2J + 1)(2J' + 1)(2F + 1)(2F' + 1)] 

V J' S \ j J' F' I \ ( F 1 F' 
J L I F J 1/1/ g -/' 



(10) 



where {} is the 6-J symbol, () is the 3-J symbol with polarization number q, and there is an 
implied normalization constant such that the cycling transition strength is 1. 

Similarly, we can calculate «2, the probability per emitted photon of leakage from the 
bright state into the dark state. Note that there are two leakage paths out of the bright 
state: via coupling to the P3/2I/ + 1/2,/+ 1/2) state (due to 7r-polarized laser light) or to 
the /3/2I/+ 1/2,/— 1/2) state (due to er~-polarized radiation). These two paths yield a 
leakage probability: 

c* = (i + - + ( 2 A 2 ) (+) 2 MaA+M2 -f- (id 

where A2 = ujkfp is the hyperfme splitting of the P3/2 levels, (P_) is the fraction of 7r 
((T _ )-polarized laser power, and (M2-) is the dipole branching ratio for 7r (ct - )-polarized 
light impurity. These branching ratios are given by: 

M 27r = C(J + 1/2, /+ 1/2; I + 1/2, 1 + 1/2) xC(J- 1/2, 1 +1/2; I -1/2, 1 +1/2) 

4/ , , 

(12) 



18/ 



and 



M 2 _ = C(I + 1/2,1+ 1/2)1 + 1/2,1 — 1/2) X C(I - 1/2,1+ 1/2; I— 1/2,1 - 1/2) 

- 167 (13) 
" 9(1 + 21)8' [L) 

2.4-2 1=1/2 case using P1/2 

For the special case of detecting the qubit state via the P\ « manifold, we need to make some 
small adjustments to our calculations. Assuming that the laser power is split roughly equally 
between the 3 polarizations, the mean number of collected, scattered photons when starting 
in one of the bright, F = 1 states is: 

. V 

A = TDT) " T n2 (14) 

25 

V 
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where 7' is the radiative linewidth of the P1/2 state. 

For the dark state leakage per emitted photon, a[, we simply have a new relevant detuning 
A[ and branching ratio Mj : 

a ' 1 -M' 1 (l + s+ (^) 2 )(^ [ ) 2 . (15) 

where the relevant detuning is now A^ = whfs + w hfpi namely the S'1/2 state hyperfine 
splitting plus the P1/2 state hyperfine splitting. For the 1=1/2 ions, all allowed dipole 
transitions have relative strength of 1/3, so the dark state leakage branching ratio is: = 
1/3 x (1/3+1/3) = 2/9. 

Also, the bright state can leak into the dark state via off-resonant coupling to one of the 
Pi/2) F = 1 states with probability per emitted photon of: 

4-M' 2 (l + s+ (^) 2 )(^) 2 (16) 

where A' 2 = ^hfp> the hyperfine splitting of the P\/2 state, and M' 2 = (1/3+1/3) x 1/3 = 2/9. 

2.5 Example calculation using 111 Cd + 

We now compare the above model to experimental data for m Cd + qubits using a PMT for 
photon detection. For this calculation and experiment we use the general detection technique 
that couples to the P 3 / 2 levels. 

In the experiment, individual 111 Cd + ions (J = 1/2) are confined in a linear RF Paul trap 
with controllable axial frequency cu z /2tt = 0.5 — 2.5 MHz 32J- The qubit is stored in the first- 
order magnetic field-insensitive "clock" states: |0, 0) = |0) and |1,0) = |1). In m Cd + , for 
cr + -polarized laser light the large hyperfine splitting of the P 3 / 2 states ensures that population 
in the 'S'i /2 1 1 , 0} state is optically pumped to the S1/2II1 1) with near unit probability so the 
P3/2 detection theory given in sec. 12.4. 11 will apply well to this case (corrections are given 
in sec. I2.6.I|) . For 111 Cd + , the relevant energy splittings are: r y/2n = 60 MHz, Ai/2n = 
13.7 GHz, A2/2-7T = 800 MHz [23 03] and we choose the laser wavelength such that the 
detuning 5 ss (figure|2J). The branching ratios are calculated from eqns.ElEl and 1 131 Mi = 
2/9 and = M2- = 1/9. Because = M2-, the it- and cr~-polarized components can 
be treated together and we can define an overall laser power impurity Pimpure = Pn + P— ■ 
Experimentally, we know the detection time, to, and we can leave 77, Pimpure, and s as free 
parameters to be determined by the fit to experimental data (although in principle one could 
independently measure these parameters). Figure^shows the good agreement between theory 
and experiment. 

2.6 Theoretical detection fidelity for 111 Cd + 

2.6.1 Detection with P3/2 levels 

To determine the theoretical limit of detection fidelity for 111 Cd + using the P 3 / 2 levels we 
choose the most ideal conditions: small laser detuning from resonance (6 — > 0) and perfect 
detection beam polarization so the bright state histogram is a true Poissonian (Pimpure = 
which means 0:2 = 0). Moreover, with perfect polarization we can lower the laser intensity 
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F=1 




Fig. 3. Relevant energy level spacings and relative electric dipole strengths for selected transitions 
in ln Cd+ 33 34 . Note that the Zeeman splitting (< 10 MHz) is small compared to the hyperfine 
splittings. 
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Fig. 4. Detection histograms using a PMT and the P3/2 detection scheme for a single 111 Cd+ ion 
prepared in: a) |0,0) = |0) dark state; and b) 1 1 , 0} = |1) bright state. For the data shown, each 
state is prepared and then measured 20,000 times. Fit is to theory from text (eqn. Inland eon. l6l 
respectively) with parameters: Tr> = 150/js, r/ = 1.4 X 10~ 3 , Pimpurc = 1.5 X 10 — 3 , and s = 0.25. 
Note that the n = 1 bin in a) includes background laser light scatter that the model does not 
include. 
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(s — > 0) to eliminate power-broadening while increasing td to maintain a sufficiently bright 
"bright" state. 

Using a discriminator at photon level d, detection fidelity of the bright state is the proba- 
bility that the ion scatters more than d photons while for the dark state it is the probability 
that it scatters d or fewer photons. The overall qubit detection efficiency is the lower of these 
two numbers and it is maximized for the two fidelities being equal. Thus, the optimal qubit 
detection fidelity is: 

d d 

F = y^Pdark(n) = 1 - y^Pbright(n) (17) 

n=0 n=0 

with Ao = sTfjrj^ and ct\ = M^^-) 2 given our assumptions mentioned above. 

To calculate the detection fidelity, F, we find a\ = 1.1 x 10~ 6 for 111 Cd + . Then we can 
choose the optimal Ao for discrimination by controlling the effective light level on the ion 
(adjusting the product std). For our current value of ?y ~ 0.001 the optimum light level 
yields Ao = 5.6 and, using the discriminator level d = 0, we obtain a qubit detection fidelity 
of F = 99.5%. In practice the optimal light level is slightly higher to aid in discriminating 
the background laser scatter from the bright qubit fluorescence. 

We can obtain a useful approximate analytic result for the theoretical qubit fidelity as 
a function of detector collection efficiency by assuming a discriminator level d = 0. The 
optimal situation is when the fidelity of bright state detection and dark state detection are 
equal: F — e~ Q i A "/»; = 1 — e~ A ° (eqn. 1170 . The term aiXo/r] is typically small so we Taylor 
expand and then take the natural logarithm of both sides to find: 

Ao+lnAo^-lnK/r?) (18) 

We wish to obtain a closed form solution for Ao as a function of ct\ and r\ since ol\ depends 
only on atomic parameters. Therefore, for the typical case of optimal Ao ~ 5 — 15 it is 
reasonable to further assume that InAo <C Ao so that the ideal light level is approximately 
Ao ~ — \n(ai/rj) — ln(ry/ai) and the approximate fidelity is: 

F approx «l-^ln(^-). (19) 

This result gives an accurate scaling of the fidelity as a function of detector collection efficiency 
(figure EJ). 

In an ideal case, we could use a cavity to surround the ion |rS5| to increase the solid angle 
of light collected and a detector with very high quantum efficiency. This might produce a 
photon collection efficiency as high as ~ 30% which, using eqn. El would yield a fidelity: 
-Fidcai = 99.997%. 

At this point we hit a fundamental limitation of fluorescence state detection when using the 
P 3 / 2 levels and a qubit stored in the clock states: even with perfect polarization the |1, 0) = |1) 
state can off-resonantly couple to the |0, 0) = |0) dark state before ever reaching the cycling 
transition. This coupling is unaffected by any external parameters and is a property of the 
atomic structure only. Given this inherent error, the maximum fidelity of direct clock state 
qubit detection is given by: 

F_ = l-P M ^ dark = l-^(^) 2 (20) 
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Fig. 5. Theoretical prediction of detection infidelity (1-F) using the P3/2 detection scheme for 
ln Cd+ as a function of total detector collection efficiency (t;). Filled circles are numerically 
calculated using can, 1171 with iterations to find the optimal light level & detection time for each 
r\\ solid curve uses the approximate analytic result of eqn. 1191 



where the factor of 4/9 comes from the squared Clebsch- Gordon coefficients for transitions 
out of the |1) state in Cd + . For ions where the -P3/2 hyperfine splitting is not much larger 
than the linewidth, the "bright" clock state is not actually very bright because it couples to 
the dark manifold so easily. Direct fluorescence state detection of clock states is basically 
impractical for these ions, but there are a number of schemes to improve detection fidelity 
including using Raman transitions to shelve one of the qubit states in a more off-resonant state 
(with the detection fidelity now possibly limited by the Raman 7r-pulse fidelity) and using 
large magnetic fields to make lst-order insensitive states with rap ^ [^Sj. For 111 Cd + , 
ujhfp ^S> 7/2, so direct clock state detection yields F max — 99.94%. Reaching this limit 
requires a high quantum efficiency detector (see eon. 119(1 . a subject we address in sec.0 

2.6.2 Detection with P\/ 2 levels 

The theoretical limit on detection fidelity when coupling to the P\/ 2 levels is the leakage from 
one of the bright states into the dark state. With no laser detuning (6 — * 0) and low light level 



(s — > 0), the relevant parameters become: X' = stdV' 



and a'n 



It is always true that > A' 2 , so a' 2 > a[ which means that there is always more leakage 
from the bright state into the dark state than vice versa. Tabled provides the relevant energy 
splittings for m Cd + , 171 Yb + , and 199 Hg + and the calculated detection fidelity for given 
values of 77, the detector collection efficiency? Note that for m Cd + the advantage of using 
the P\/ 2 detection scheme instead of the P 3 / 2 scheme is only realized for large values of 77. 



3 Individual Ion Detection Using a CCD 

6 Note that the fidelity for detecting the qubit state in 171 Yb + is slightly more complicated because there is an 
allowed decay from the excited P1/2 state to a low-lying -D3/2 state which can be re-pumped to the S1/2 level 
via the [3/2] 1 / 2 level 1371 . This repumping step is identical to the standard detection step, but its infidelity 
can be ignored beacuse the off-resonant coupling is much smaller than for the Pi/ 2 detection step due to the 
energy splittings: 7[3/2]/2tt = 9.5 MHz, A H fd/2tt = 0.86 GHz, A HP [ 3 /2]/27r = 2.5 GHz, and A = 935.2 nm. 
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3.1 CCD technical overview 

To benefit from the spatial resolution of a CCD one must typically use an intensified CCD to 
obtain a signal that is much larger than the CCD readout noise (eqn. I21fl f In an intensified 
CCD, single photons incident on the front-screen photocathode produce electrons that are 
accelerated across a multichannel plate to start a localized electron avalanche that impacts 
a phosphor screen (figure EJ) . Negative (positive) biasing of the photocathode relative to the 
multichannel plate produces significant (no) electron acceleration; this differential effect allows 
for rapid gating of the intensifier and reduces background counts. Visible wavelength photons 
emitted from the phosphor are then coupled (via lenses or fiber optics) to a standard CCD 
that converts the photons to charge for readout. The charge on each pixel of the CCD is then 
measured with the value proportional to the incident light intensity. A computer PCI board 
converts the analog voltage signal for each pixel to a digital integer for computer processing. 




Fig. 6. Schematic diagram of an intensified CCD camera imaging tube. Incident photons impact 
the photocathode where they are converted to electrons and accelerated across the multichannel 
plate. These accelerated electrons strike the phosphor screen and produce visible photons that are 
coupled into fiber optics and fed onto a standard CCD for readout. Diagram used with permission 
of Princeton Instruments / Acton. 

Technical noise considerations dictate the important CCD characteristic of readout speed. 
For a qubit in the bright state, the average number of electrons produced at the photocathode 
by photons incident on the detector during detection is: An = rii on r]r)Tdfl/ATr where rii on is 
the number of photons emitted by the ion, rju is the quantum efficiency of photon to electron 

c Recent advances in electron-multiplying CCD (EMCCD) architectures allow for single-photon detection at 
high readout speeds with low noise (although with no gating capability) 1401 . The following discussion, in 
particular eqn. 1211 is equally valid for EMCCDs by using the appropriate gain and quantum efficiency. 

Table 1. Relevant energy splitting parameters for the / = 1/2, P\/2 detection scheme using 
ln Cd+ 171 Yb+ and 199 Hg+ with detection fidelities calculated using cqn.^for 

varying detector collection efficiency r). 





ln Cd+ 


171 Yb + 


199 Hg+ 


A(S 1/2 - P 1/2 ) (nm) 


226.5 


369.5 


194 


i /2tt (MHz) 


50 


23 


70 


Ahfs/2tt (GHz) 


14.5 


12.6 


40.5 


A^ Fp /27r (GHz) 


2 


2.1 


6.9 


n = 0.001 


F = 96.7% 


F = 99.33% 


F = 99.43% 


n = 0.01 


F = 99.65% 


F = 99.93% 


F = 99.943% 


n = 0.3 


F = 99.988% 


F = 99.998% 


F = 99.998% 
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conversion, T is the optical transmission between the ion and the detector, and ^ is the 
solid angle of light collection. The remaining stages of electron intensification, conversion to 
photons, and then photon-induced charge production can be summarized by a gain factor 
g so that the CCD rms shot noise will be: g\f\o- To readout the charge on a pixel, the 
CCD controller must first clear the charge accumulated during the previous pixel readout. 
The imperfect repeatability of this process will induce some noise on the signal, with faster 
readout speeds leading to less perfect charge cleaning and therefore more noise. This rms noise 
per readout, r, is uncorrelated to the shot noise so the two noise sources add in quadrature to 
produce a total noise: ^/g 2 Ao + (fcr) 2 where fc is the total number of pixels readout. Therefore, 
the total signal-to-noise ratio (SNR) becomes: 



Many CCDs allow on-chip binning of multiple pixel charges together before readout. This 
on-chip binning increases readout speed and decreases readout noise (because less pixels are 
readout) at the expense of decreasing spatial resolution (although in principle resolutions 
above one pixel per ion are unnecessary for distinguishing between the bright and dark states). 
If VAo <C kr/g then the signal is readout noise limited and SNR w so the SNR will 
improve linearly with the number of pixels binned. For timescales typical of an ion fluorescence 
experiment, fcr « 10 and so using an intensifier with g ^> 1 allows one to use the gain to 
overwhelm the readout noise. In this shot noise limited regime, the signal-to-noise ratio 
becomes SNR w \/Ao- 

3.2 CCD experimental usage 

For efficient experimental detection we use a 28 x 28 pixel box on-chip binned 4x4 such that 
there is an effective box size of 7 x 7 "super" -pixels. This box size was chosen to match the 
imager-magnificd ion-ion spacing in a trap with uj z /2tt = 2.0 MHz so that each ion would 
be centered in its box with no overlap of boxes; the 4x4 binning improves readout speed 
yet still offers enough imaging resolution to provide real-time monitoring of ion/imaging 
system drifts. The detection signal is formed from the integrated electron counts of these 
49 pixels minus the constant offset due to the non-zero readout charge maintained on each 
CCD pixel. Note that this offset is not the same as the negligible CCD dark counts that are 
caused by inadvertent electron transfer to the back-screen due to thermal effects. We are able 
to significantly suppress this thermal effect by cooling the CCD. When histogrammed, the 
integrated counts for n incident photons shows a distribution because multiple electrons and 
photons are produced for each incident photon with a mean number of integrated counts for 
this box size of ~ 100 integrated counts per incident photon. Experimentally, the optimal 
discrimination level is chosen by equalizing the fraction of misidentified states, i.e. choosing 
a discrimination level such that an equal fraction of 1 1,0) = |1) is misidentified as "dark" 
and 1 0,0) = |0) is mislabeled as "bright". Using this optimal discriminator the efficiency of 
detecting a single ion's quantum state is > 99.4% with the lower detection fidelity attributable 
to the CCD's smaller light collection angle (the PMT counts all the incident light while the 
CCD only uses counts inside the ion pixel box). 



SNR = 




(21) 
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4 Multiple Ion Detection Using a CCD 
4-1 Experimental implementation 

Of more practical interest than single-ion detection is the ability to detect multiple ion quan- 
tum states simultaneously using the CCD's spatial resolution. For each ion we select a box 
or region of interest (ROI) that determines the pixels over which the CCD will integrate to 
determine the ion's state. Ion ROI box edges are positioned on adjacent pixels with no overlap 
to maximize detection fidelity and on-chip binning is used to enhance readout speed. 

Figure [7^ and b shows detection histograms and the actual CCD "pictures" for three ions 
prepared in the |000) dark state and the | 111) bright state. In fig.0: we show the detection 
histograms and post-selected examples of all possible pictures for an equal superposition 
state achieved by preparing the state |000) and then applying a microwave 7r/2-pulse at 
the hyperfine splitting (x>hfs/27t = 14.5 GHz equally on the three ions to produce the state 
(|0> + |1»® 3 . 

4-2 Possible additional errors 

In addition to the single-ion detection errors due to off-resonant coupling (secQ, simultaneous 
multiple ion detection on a CCD can produce some additional errors. One possible error is 
that the detection beam has a waist of ~ 10 fim, compared to an ion-ion spacing of ~ 4 /zm 
for this experiment, leading to unequal illumination of the ions (visible in the histograms of 
fig. where the middle ion is brighter than the outer two ions). In general, this unequal 
illumination results in different constant offsets for each ion's light level that are easy to 
correct. 

The more significant error source is the potential for adjacent ion effects, namely opti- 
cal pumping from ion-to-ion and box-to-box leakage on the CCD. Pumping an ion into the 
"bright" state from the fluorescence of the adjacent state is completely neglible compared to 
the effect of the laser. The light intensity due an adjacent bright ion is: h ^J.'i where u> is 
the frequency of the radiated light and x is the inter-ion spacing. In comparison, the laser 
intensity is typically: /i aS cr ~ -4 at = where A is the wavelength of the radiated light. 

Thus, the light intensity effect of the adjacent ion relative to the laser is: 

= 4^- < 22 > 

For 111 Cd + , A = 214.5 nm and with three ions in a trap with ui z /2ir w 850 kHz, x w 4 /xm. 
Thus, iion/iianer ~ 7 x 10 -4 , which is completely negligible. 

Box-to-box leakage on the CCD can be most clearly seen by examining the conditional 
probabilities of the ions dependent on the state of the other ions. By independently rotating 
the qubits with microwave pulses, there should be no correlation between qubit states, but 
light leakage between boxes would produce such correlations. Using the data of figure [7J 
we find that adjacent ion pairs (left/center and center/right) exhibit conditional probability 
correlations of ~ 1.2%. With this leakage, the detection fidelity for individual qubits in the 
presence of other qubits is ~ 98%. This fidelity is fully consistent with a small amount of 
leakage of the light from one ion into the ROI of the adjacent ion, an effect most likely caused 
by ion/imaging system drifts over a few minutes. Further mechanical stabilization of the 
imaging system should eliminate this leakage. 
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Dark state data 



Bright state data 
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Fig. 7. Detection histograms for three ion qubits in the states: a) 1 0, 0) = |0); b) |1,0) = |1); and 
c) an equal superposition state prepared by starting in |000) and applying a microwave 7r/2-pulse 
equally on all three qubits to produce the state (|0) + |1))® 3 . Each graph contains 4000 trials. 
Note that ion #2 has the most integrated counts due to unequal illumination by the detection 
beam. Adjacent to each histogram are examples of the post-selected single-shot images acquired 
by the CCD for each case: a) all "dark" (|000)); b) all "bright" (|111)); and c) all combinations 
of "dark" and "bright" . The dashed white lines indicate the boundaries of the regions of interest 
used to determine the qubit state via integrated CCD counts. 
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5 Future Work 

The current detection fidelity of ~ 98% is highly efficient, but improvements in light collection 
and/or CCD quantum efficiency would be necessary to increase the fidelity to perform fault- 
tolerant quantum computing with a reasonable number of qubits. In addition, the CCD total 
readout time of ~ 15 ms is much longer than a typical gate time of r gate ~ 100/is, limiting 
the CCD's effectiveness for algorithms with steps conditional upon state detection (such as 
quantum error-correction |19| or the quantum Fourier transform |17|). Both limitations are 
ultimately technical. Current state-of-the-art unclassified CCDs operating in the near-IR 
can produce readout speeds of ~ 10 MHz per pixel or ~ 2.5 /is per ion ROI with quantum 
efficiencies of ~ 60%. With these types of performance characteristics, employing CCDs for 
simultaneous multiple ion detection in feedback algorithms will be practical and efficient. 
Note that the fundamental limit on detection speed is given by the lifetime of the excited 
state: r = — . Thus, with high-efficiency light collection, the state detection time for Cd + 
could be as short as ~ 50 — 200 ns. 

Large-scale ion trap arrays have been proposed [5] an d small, scalable traps have 
been successfully produced using microfabrication techniques 02 S3] ■ Integrating multiple 
qubit detection via the CCD with these traps would produce a truly scalable qubit processing 
architecture. 
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